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Abstract:  This  paper  presents  a  novel  method  for  the  element  characterization  of  printed  re¬ 
flectarray  antennas.  Their  characterization  requires  the  calculation  of  the  phase  of  the  reflection 
coefficient  versus  the  patch  dimension  (phase  diagram),  when  the  structure  is  illuminated  by  a 
uniform  plane  wave.  The  calculation  of  the  phase  diagram  is  based  on  the  full-wave  analysis  of  the 
structure  in  a  limited  number  of  points,  and  a  polynomial  interpolation  of  the  phase,  which  takes 
advantage  of  a  perturbational  technique  for  estimating  the  phase  derivative.  This  approach  sig¬ 
nificantly  reduces  the  computing  time  for  the  calculation  of  the  phase  diagram,  as  shown  through 
an  example. 
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1  Introduction 

Printed  reflectarrays  antennas  received  a  lot  of  attention  in  the  last  decade  as  a  valid  alternative 
to  classical  reflector  antennas,  for  both  scientific  and  commercial  applications  [1]— [5] .  Similarly  to 
reflector  antennas,  reflectarrays  comprise  a  low-gain  feed,  located  in  front  of  a  reflecting  surface. 
The  peculiarity  of  reflectarrays  is  that  the  reflector  is  a  planar  array  of  printed  patches  on  a 
grounded  dielectric  substrate  (Fig.  1).  In  order  to  obtain  the  focussing  effect,  the  patches  have 
slightly  different  size  and/or  shape,  thus  leading  to  a  proper  variation  of  the  phase  of  the  local 
reflection  coefficient  [1]. 

Different  techniques  have  been  proposed  for  phase  shifting  the  local  reflection  coefficient: 
among  them,  the  most  popular  are  based  on  the  use  of  identical  microstrip  patches  with  dif¬ 
ferent  length  phase-delay  lines  [1,  3],  variable  size  microstrip  patches  [2,  5],  and  identical  patches 


reference  plane 


Figure  1:  Principle  of  operation  of  a  printed  reflectarray. 
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Figure  2:  Classical  shapes  used  in  printed  reflect  arrays. 


with  different  rotation  angles  [4]  (Fig.  2).  The  choice  of  the  element  shape  is  a  key  issue  in  the 
design  of  reflectarray  antennas.  The  performance  of  different  element  shapes  have  been  compared 
in  a  systematic  way  in  [6]. 

The  major  advantages  of  printed  reflectarray  antennas  are  low  profile,  reduced  weight,  and  easy 
manufacturing.  Moreover,  especially  in  spaceborne  systems,  printed  reflectarray  can  be  folded  for 
launching  and  deployed  for  operation.  The  only  significant  disadvantage  is  the  narrow  bandwidth 
of  these  antennas,  typically  limited  to  few  percent. 

In  the  design  of  printed  reflect  arrays,  geometrical  optics  is  usually  applied  for  the  calculation 
of  the  electrical  length  of  the  ray  path  from  the  feed  to  the  reflector  and  from  the  reflector 
to  the  reference  plane  (Fig.  1).  Conversely,  a  full-wave  analysis  is  applied  for  the  calculation 
of  the  phase  shift  due  to  the  reflection.  This  analysis  is  typically  performed  under  the  local 
periodicity  approximation:  each  metal  patch  is  considered  in  a  periodic  array  of  identical  elements, 
thus  neglecting  the  variation  of  mutual  coupling  due  to  the  small  difference  between  neighboring 
elements.  The  output  of  the  full-wave  analysis  is  the  “phase  diagram”,  which  is  the  plot  of  the 
phase  of  the  reflection  coefficient  versus  the  geometrical  dimension  A  of  the  array  element  (Fig.  3). 
Once  the  path  lengths  and  the  phase  diagram  have  been  calculated,  the  proper  element  dimension 
is  selected,  in  order  to  create  the  required  radiation  pattern  of  the  reflectarray  antenna. 

Whichever  is  the  adopted  full-wave  approach,  the  calculation  of  the  phase  diagram  is  the 
most  CPU  intensive  part  of  the  analysis.  In  this  paper,  we  propose  the  use  of  an  interpolation 
technique,  which  is  based  on  a  perturbational  approach  in  conjunction  with  an  integral  method, 
named  the  MoM/BI-RME  method.  This  technique  allows  for  a  substantial  reduction  in  CPU 
time. 
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Figure  3:  Example  of  a  phase  diagram:  a)  geometry  of  the  patch  and  geometrical  parameter  A; 
b)  phase  of  the  reflection  coefficient  versus  A  ,  in  the  case  of  normal  incidence,  vertical  polarization. 


April  19-23,  2004  -  Syracuse,  NY  ©  2004  ACES 


20th  Annual  Review  of  Progress  in  Applied  Computational  Electromagnetics 


2  Phase  diagram  calculation  by  the  MoM/BI— RME  method 

In  our  approadi,  the  determination  of  the  phase  diagram  is  based  on  the  MoM/BI-RME  method. 
This  method  consists  in  the  solution  of  an  electric-held  integral  equation  by  the  Method  of  Mo¬ 
ments  (MoM)  with  entire-domain  basis  functions.  The  vector  basis  functions  are  related  to  scalar 
functions  which  are  the  first  N  eigenfunctions  of  the  Helmholtz  equation  in  a  domain  S  coin¬ 
cident  with  the  cross-section  of  the  patch 

V2^  +  k2  if;  =  0  in  S'  (1) 

with  Dirichlet  =  0)  or  Neumann  (dyj jdn  =  0)  boundary  condition  on  dS  (n  is  the  outward 
normal  unit  vector  on  dS).  The  eigenfunctions  ? bi  and  the  corresponding  eigenvalues  Ki  are 
numerically  calculated  by  the  Boundary  Integral-Resonant  Mode  Expansion  (BI-RME)  method, 
which  applies  to  printed  elements  with  an  arbitrary  shape  [7]. 

The  MoM/BI-RME  method  was  originally  developed  for  the  analysis  of  frequency  selective 
surfaces  (FSS)  [8,  9].  In  that  case,  the  method  provides  the  frequency  response  of  the  FSS  for  a 
given  geometry.  Its  application  requires  two  different  steps:  the  determination  of  the  basis  func¬ 
tions,  which  is  frequency  independent,  and  the  solution  of  the  MoM  problem,  which  is  performed 
frequency-by-frequency  (Fig.  4a).  The  major  advantage  of  this  method  derives  from  the  use  of 
entire-domain  basis  functions,  which  must  be  calculated  only  once  and  lead  to  a  MoM  matrix 
problem  with  a  small  number  (IV)  of  unknowns. 


a  b 

Figure  4:  Flow-chart  of  the  MoM/BI-RME  method:  a)  calculation  of  the  frequency  response  of 
frequency  selective  surfaces;  b )  calculation  of  the  phase  diagram  of  printed  reflectarrays. 
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The  typical  CPU  time  for  the  calculation  of  the  frequency  response  of  FSS  in  100  points  is  in  the 
order  of  5  4-  6  seconds.  This  CPU  time  is  almost  equally  balanced  between  the  calculation  of  the 
basis  functions  and  the  frequency  loop  [8]. 

The  analysis  of  reflectarrays  requires  a  slightly  different  approach.  In  fact,  the  calculation  of 
the  phase  diagram  is  performed  at  a  single  frequency,  by  varying  the  geometry  of  the  patch.  This 
implies  that  the  entire-domain  basis  functions  must  be  calculated  for  each  value  of  the  geometrical 
parameter  A  (Fig.  46),  thus  weakening  the  advantage  of  using  the  MoM/BI-RME  method. 

An  example  of  phase  diagram  calculated  by  the  MoM/BI-RME  method  is  reported  in  Fig.  3, 
in  the  case  of  a  patch  element  with  ridged  rectangular  shape. 


3  Interpolation  method  and  perturbational  technique 

A  significant  reduction  in  computing  time  can  achieved  by  using  an  interpolation  technique.  Since 
the  phase  diagram  is  typically  a  smooth  function  of  the  geometrical  dimension  A  (see  Fig.  3),  it 
could  be  well  fitted  by  using  interpolating  polynomial  functions.  In  particular,  the  range  of 
variation  of  A  is  divided  into  a  number  of  identical  intervals.  In  each  interval  [A,  A],  the  phase 
diagram  is  represented  as  a  third  order  polynomial  function 

0(A)  =  cq  +  Ci  A  +  C2  A"  +  C3  A3  in  [A,  A  ]  (2) 


and  the  coefficients  Co,  . . . ,  C3  are  determined  from  the  knowledge  of  the  phase  0  and  its  derivative 
d(j)/dA  in  A  and  A.  It  is  found  that  few  intervals  (usually  less  than  ten)  are  needed  for  an 
accurate  interpolation  (error  smaller  than  1  deg).  The  use  of  this  technique  significantly  reduces 
the  number  of  analyses,  but  requires  an  additional  effort  for  the  calculation  of  the  phase  derivative. 
The  derivative  could  be  obtained  numerically  by  the  calculation  of  the  phase  in  two  close  points, 
namely  A  and  A'  =  A(1  +  e),  with  |e|  <C  1.  Nevertheless,  even  in  this  case  the  use  of  interpolating 
functions  is  advantageous  over  the  standard  approach. 

Another  significant  improvement  in  computational  efficiency  can  be  accomplished  by  adopting 
a  perturbational  technique  in  conjunction  with  the  MoM/BI-RME  method.  In  fact,  being  0j  the 
functions  on  a  patch  with  cross-section  S  (corresponding  to  a  given  value  of  A),  functions  0)  in 
S'  (corresponding  to  A')  can  be  obtained  as  a  linear  combination  of  the  N  functions  0* 

N 

0'  =  ^  Xij  0*.  (3) 

i— 1 


The  unknown  coefficients  xt]  are  obtained  as  follows.  A  variational  formula  for  k'  (which  represent 
the  eigenvalues  of  the  Helmholtz  equation  in  the  domain  S')  is  used  [10] 


K 


/  2 


IS' 


V0'  •  V0'  dS 


di 


(4) 


where  (3  =  1  or  (3  =  0  in  the  case  of  Dirichlet  or  Neumann  boundary  condition,  respectively.  By 
substituting  (3)  into  (4),  and  applying  Rayleigh- Ritz  method  to  find  the  stationary  value  of  the 
resulting  expression,  a  generalized  linear  eigenvalue  problem  with  dimension  N  is  obtained.  The 
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Figure  5:  Analysis  of  the  ridge-shape  element  by  the  perturbational  technique:  a)  geometry  of  the 
element  and  perturbed  region;  b)  phase  diagram  of  the  ridge  element,  calculated  by  the  standard 
method  (dots)  and  by  the  interpolation  procedure  (solid  line). 


solution  of  this  problem  yields  the  set  of  n'-  as  eigenvalues,  and  coefficients  xtj  as  the  corresponding 
eigenvectors.  Since  the  dimension  N  of  this  eigenvalue  problem  is  small,  the  new  basis  functions  -0' 
can  be  obtained  with  a  limited  computational  effort,  much  smaller  than  the  standard  calculation 
by  the  BI-RME  method. 

The  perturbational  approach  can  be  also  exploited  in  the  MoM  portion  of  the  code  (Fig.  46). 
In  fact,  the  coupling  integrals  in  the  modified  domain  S'  (between  the  basis  functions  and  the 
Floquet  modes)  can  be  written  as  a  linear  combination  of  the  coupling  integrals  in  the  original 
domain  S,  plus  line  integrals  on  the  modified  portion  of  the  boundary.  In  such  a  way,  the  most 
CPU-intensive  part  of  the  application  of  the  MoM  can  be  completely  skipped,  and  only  the 
solution  of  a  (small)  linear  system  must  be  performed. 

4  Numerical  results 

As  an  example,  the  phase  diagram  of  the  ridge  element,  already  determined  by  the  standard 
technique  (Fig.  36),  is  calculated  by  using  the  MoM/BI-RME  method  in  conjunction  with  the 
perturbational  approach.  Fig.  5a  shows  the  structure  and  the  portion  of  the  boundary  which  is 
displaced  for  the  calculation  of  the  phase  derivative.  The  phase  diagram  is  reported  in  Fig.  56. 
The  analysis  was  performed  in  only  four  points  (squared  markers),  where  the  reflection  phase  is 
calculated  by  the  MoM/BI-RME  method  and  the  phase  derivative  is  obtained  by  the  perturba¬ 
tional  technique.  The  interpolated  phase  diagram  (solid  line)  practically  coincides  with  the  one 
calculated  in  49  points  by  the  standard  method  (dots).  The  CPU  time  was  18  sec  by  the  new 
approach,  instead  of  142  sec  by  the  standard  method  (on  a  PC  Pentium  IV  @  1.4  GHz). 
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5  Conclusion 

This  paper  presented  a  novel  technique  for  the  fast  modelling  of  printed  reflectarray  elements 
with  arbitrary  shapes.  The  proposed  approach  combines  the  MoM/BI-RME  method  (used  for 
the  full-wave  analysis  of  the  structure)  with  a  perturbational  technique,  which  allows  for  obtaining 
an  interpolated  plot  of  the  phase  diagram  in  a  significantly  reduced  computing  time.  An  example 
demonstrated  the  effectiveness  of  the  proposed  method. 
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